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It has been demonstrated experimentally that non-axially symmetric vortices precess around the 
centre of a Bose-Einstein condensate. Two types of single vortex states have been observed, usually 
referred to as the S-vortex and the U- vortex. We study theoretically the single vortex excitations in 
spherical and elongated condensates as a function of the interaction strength. We solve numerically 
the Gross-Pitaevskii equation and calculate the angular momentum as a function of precession 
frequency. The existence of two types of vortices means that we have two different precession 
frequencies for each angular momentum value. As the interaction strength increases the vortex 
lines bend and the precession frequencies shift to lower values. We establish that for given angular 
momentum the S-vortex has higher energy than the U-vortex in a rotating elongated condensate. 
We show that the S-vortex is related to the solitonic vortex which is a nonlinear excitation in the 
nonrotating system. For small interaction strengths the S-vortex is related to the dark soliton. In 
the dilute limit a lowest Landau level calculation provides an analytic description of these vortex 
modes in terms of the harmonic oscillator states. 



I. INTRODUCTION 

Quantised vortices in superfluids have been the sub- 
ject of an important research area since their theoretical 
prediction in the context of liquid helium p] . In the last 
years there has been intense experimental and theoretical 
study on the creation and dynamics of single vortices as 
well as vortex lattices in ultracold atomic Bose-Einstein 
condensates (BECs) 0,0. These systems offer the pos- 
sibility for direct observation of vortex lines and of their 
dymamics. They are thus excellent systems for a com- 
prehensive theoretical study of vortex lines in three di- 
mensions. 

In an experiment performed in a spherical trap 0| sin- 
gle vortices were produced by a technique which induces 
a 2tt phase winding using a rotating laser beam. The 
vortices were located off the trap center, and they were 
seen to precess robustly around it at an almost constant 
angular velocity. Surprisingly the precession frequency 
was found to be almost independent of the displacement 
of the vortex from the trap center. 

In another important experiment it was established 
that the vortex lines in the finite trap geometry can have 
nontrivial shapes ■ It was observed that a single vortex 
line may have an S shape or a U shape. The length and 
precise shape of the vortex can vary and this is directly 
related to the angular momentum it carries. Theoretical 
work in high density very elongated condensates similar 
to those of Ref. 0] showed that the U-shaped vortex can 
be the ground state of the rotating system |(| . Some 
differences between the S and the U-vortex were reported 
in0. 

In this paper, we establish that an axisymmetric BEC 
generically supports both S and U-shaped vortices. We 
study these single vortices as precessing stationary states 



within the nonlinear Gross-Pitaevskii theory for a range 
of interaction strengths. Their frequency of precession 
around the trap center depends critically on the features 
of the vortex line. Precession frequencies for U-vortices 
depend only weakly on their angular momentum and fall 
in a relatively narrow frequency range which depends on 
the interaction strength. S-vortices have higher energy 
and their precession frequencies depend strongly on their 
angular momentum. The nonrotating member of this 
family is a non-axially symmetric vortex solution of the 
nonlinear theory. This has been called a solitonic vor- 
tex; it breaks spontaneously the symmetry in an axially 
symmetric potential and pierces through the trap in a 
direction perpendicular to the symmetry axis [8|, |9j . 

In the dilute limit the low-energy states can be de- 
scribed within a lowest Landau level (LLL) ansatz [To| . 
Motivated by the fact that experiments are usually per- 
formed in elongated condensates we generalize the LLL 
ansatz to include longitudinal modes. Both types of sin- 
gle vortex states can then be described within this for- 
malism and some simple analytic expressions capture the 
main features of the vortex modes. Low density conden- 
sates close to the LLL regime have been created recently 
by fast rotation where a vortex lattice was formed [TlT ] . 

The outline of the paper is as follows. In Sec. II we 
formulate the problem and find numerically two types of 
single vortex solutions of the Gross-Pitaevskii equation. 
We apply our method for a spherical trap as in |4| as well 
as for elongated condensates. In Sec. Ill we formulate the 
problem in the LLL augmented to include the longitudi- 
nal modes and we study the S and U-vortices within this 
formalism. Finally, Sec. IV contains a summary and our 
concluding remarks. 



II. PRECESSING VORTICES 

Motivated by the experiment |4J where off-center vor- 
tices precess at an almost constant angular velocity 
around the center of the condensate, we look here for 
vortices which are stationary states in a frame rotating 
at a constant frequency. In this paper we suppose that 
the energy of a BEC is given by the Gross-Pitaevskii 
theory. We consider an axially symmetric harmonic po- 
tential with frequencies loj_ in the transverse direction 
and uj\\ in the longitudinal direction and implement ra- 
tionalised units by the substitutions 



a± x, 



a±y, 



a \\ z l 



N 1 ' 2 



1/2 



where "J is the condensate wave function, N is the num- 
ber of atoms in the condensate, and dj_ 11 = y/h/mu>±^ 
are the oscillator lengths in the transverse and longitu- 
dinal directions. The wave function normalisation condi- 
tion reads 



\WdV = 1. 



(2) 



The energy functional per particle in units of tku±_ has 
the form 
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where p, z, and <f> are cylindrical coordinates, is the 
gradient on the transverse plane, a s is the scattering 
length, and the parameter 



(4) 



is the ratio of the two frequencies. The value (3 = 1 
corresponds to a spherical trap. The quartic term gives 
the strength of the interactions and we introduce for later 
convenience the parameter 



S = 



1 Na, 



2V2tt an 



(5) 



which is a dimensionless measure of the interaction 
strength relative to the transverse trapping energy Juo_i_. 

Stationary states which precess at a constant angular 
frequency lo are extrema of the energy in the rotating 
frame: 



E IC 



E-lo£, 



(6) 



where £ is the angular momentum per particle along the 
symmetry axis, in units of h: 



(7) 
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FIG. 1: The particle density |\I/| 2 for vortex solutions of 
the Gross-Pitaevskii model in a spherical trap (/3 = 1) at 
S — 85. The parameters correspond to the conditions of 
the experiment Q. We show the plane containing the vor- 
tex line and the z axis along which the angular momentum 
i is measured. The images correspond to four representative 
values of I = 0.17, 0.49, 0.73, 0.98 with precession frequencies 
lo = 0.268, 0.221, 0.205, 0.196 respectively, in units of lo±. The 
size of the frames shown is 12x12 in units of a±. White colour 
corresponds to no particles and black colour to maximum par- 
ticle density. 



The star denotes complex conjugation, and the frequency 
lo is measured in units of lo±. 

Extrcmizing numerically the energy functional 10 is 
often not practically convenient because some of the in- 
teresting solutions may not correspond to its minima. 
Therefore, we dehne the Lyapunov functional 



Kot 



a i 



E+^(£-h) 2 , 



(8) 



where we have introduced the constants a\ , b± . An ex- 
tremum of E' mt corresponds to an extremum of E rot with 
lo = a\{b\ — £). We find the minima of E' mt by using a 
variant of a numerical norm-preserving relaxation algo- 
rithm |12| which, in effect, capitalizes on a virial relation 
in order to converge to a wave function of unit norm. 
The method is applied on a three-dimensional grid in 
Cartesian coordinates. The advantage of using E' mt is 
that we can find stationary points which are not neces- 
sarily minima of E rot . The constants a\,b\ should be 
chosen so as to ensure convergence to the desired value 
of the angular momentum, say £q. This usually amounts 
to choosing b\ close to £q, while a± should be positive and 
large enough so that the stationary point at £q be a mini- 
mum. Once the solutions for a range of angular momenta 
are known the frequency is given by lo — dE/d£. 

We first apply the method for parameters similar to 
those in Ref. jj]: j3 = 1 and 5 = 85. Since the trap is 
spherical any arbitrary axis in space can be chosen as 
the z-axis along which £ is measured. We find vortex 
states for the range of angular momenta < £ < 1 and 
we present images of the particle density for some repre- 
sentative ones in Fig. ^ The condensate is at the non- 
rotating spherically symmetric ground state for £ — 0. 
For £ > a vortex line enters the condensate. As £ is 
increased the vortex line is closer to the center and it is 
curved for £ < 1; it becomes the axially symmetric vortex 
at the center of the condensate for I = 1. Following 0)0 
we call the vortex mode in Fig. ^ the TJ -vortex. 

The angular momentum £ as a function of the preces- 
sion frequency lo of the vortex around the z-axis is shown 
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FIG. 2: Angular momentum as a function of the precession 
frequency for f) = 1, 6 = 85 (solid line) which corresponds to 
Fig- Q (and to experiment |4)). We have a U- vortex branch 
with precession frequencies 0.20 < ui < 0.35 in units of ui±. 
We include the corresponding line (dashed) for a smaller in- 
teraction strength 8 — 41 for comparison. 



in Fig. |21 We find that single vortex states with < I < 1 
have frequencies in the range 0.19 < u> < 0.35. However, 
except for very low values of £ the branch has frequen- 
cies in a narrow frequency range, say, uj w 0.19 to 0.25. 
In physical units the frequency / defined from lu = 2nf 
takes values in the range 1.56 Hz < / < 1.95 Hz (e.g., 
/ = 1.72 Hz for £ = 0.5); which is in agreement with 
the value 1.8 Hz measured in 0. In Ref. [l4| the pre- 
cession frequency of the vortex is explained in terms of 
an anomalous mode of the Bogoliubov spectrum and it is 
related to its stability. The frequencies 1.58 Hz given in 
the above reference should be compared with the limiting 
frequency (1.56 Hz) obtained in our calculation. A some- 
what counterintuitive feature of the results presented in 
Fig. [21 is that the precession frequency decreases as the 
angular momentum increases, which means that vortices 
closer to the axis precess slower. This point is discussed 
further at the end of this section in connection with sta- 
bility. 

Further in this section we expand our calculations to 
axially symmetric elongated traps. This is an impor- 
tant direction since the majority of experiments involv- 
ing vortices have been performed in elongated conden- 
sates. Previous theoretical work on few-vortex solutions 
has focused on very dense elongated condensates which 
are practically in the Thomas- Fermi regime 0, IE an d 
satisfy the condition 6 ^> 1, /3. In this section we shall 
focus on the regime S ~ 1. 

We present here a calculation for (3 = 1/2 and 5 = 2. 
Using Na atoms and trapping frequencies u>± = 27r x 7.8 
Hz, = ll>±/2, these values are obtained for N = 27000 
atoms (N = 6700 for Rb). In Fig.[3|we show by a solid 
line the angular momentum for the U- vortex as a function 



FIG. 3: Angular momentum as a function of the precession 
frequency for vortices in an axially symmetric elongated trap 
with (5 = 1/2, 5 = 2. The solid line shows the U- vortex 
branch. The dashed line is the S-vortex branch which contains 
the solitonic vortex for £ — (at u) = 0). Both modes give 
the axially symmetric vortex for I — 1 but this is obtained at 
different limiting frequencies: u = 0.64 for the U-vortex and 
uo = 0.44 for the S-vortex. 
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FIG. 4: The particle density |<1/| 2 for the U-vortex correspond- 
ing to the solid line in Fig. [3](/3 = 1/2, 8 = 2) for four values 
of the angular momentum £. We show here the plane con- 
taining the symmetry axis and the vortex line. At £ = we 
have the nonrotating ground state and for t = 1 we obtain 
the axially symmetric vortex. For < £ < 1 the vortex lines 
are slightly curved. The size of the frames shown is 6x10 (in 
units of a±). 



of the precession frequency. The angular momentum is 
a decreasing function of u>, and the precession frequency 
range over which the U-vortex exists is narrow and it is 
located at much higher frequencies compared to Fig. [2] 
In Fig. 0] we present images of the U- vortices for various 
values of the angular momentum. Vortices are off-center 
and slightly curved. The axially symmetric vortex is ob- 
tained as the limit for a frequency lojj = 0.64. 

Refs. HHEi study a solution of the Gross-Pitaevskii 
equation which has been called the solitonic vortex. This 
appears to be different than all the vortices discussed in 
this section so far. The solitonic vortex is an excited state 
in an elongated non-rotating condensate which pierces 
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FIG. 5: The particle density |^| 2 for the S- vortex corre- 
sponding to the dashed line in Fig. [3] (/? = 1/2, 5 = 2) for 
four values of the angular momentum I. We show here the 
plane containing the symmetry axis and the vortex line. At 
I — we have the solitonic vortex and for I = 1 we obtain 
the axially symmetric vortex. For < I < 1 the vortex lines 
are slightly curved. The size of the frames shown is 6x10 (in 
units of a±). 
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through the trap in a direction perpendicular to the sym- 
metry axis and therefore it has angular momentum £ = 
(along the symmetry axis). An example is shown in the 
first entry of Fig. The solitonic vortex exists for suffi- 
ciently high interaction strengths and it bifurcates from 
the dark soliton solution. 

We notice that the solitonic vortex breaks sponta- 
neously the axial symmetry and if we assume a finite 
precession frequency uj, it may acquire a nonvanishing 
angular momentum i = £(oj). In order to explore this 
possibility we use the solitonic vortex solution as an ini- 
tial guess in our numerical algorithm and choose appro- 
priate values for the constants in Eq. 10 ■ Representative 
results for vortex solutions with < I < 1 are shown in 
Fig- El for an elongated condensate (fi = 1/2) and for an 
interaction strength constant 5 = 2. The first entry, for 
I = 0, is the solitonic vortex and the next entries show 
almost straight vortex lines which are tilted with respect 
to the z axis. As the frequency increases the vortex line 
progressively aligns with the symmetry axis while its an- 
gular momentum increases. The dependence of £ on w is 
shown by a dashed line in Fig. This mode, which may 
be called the S-vortex following [7j, gives the axially sym- 
metric vortex (£ = 1) at a limiting frequency = 0.44 
(this value should be compared with Eq. (fH?|) in Sec. III.) 
Thus the S and the U-vortex are two different types of 
single vortex states in elongated condensates. 

The picture summarized in Figs.|3^ and [3] is appar- 
ently robust for large interaction strengths. In partic- 
ular, in the experiments in Q and in Refs. 0, IE S 
and U-vortices are found in condensates which are prac- 
tically in the Thomas- Fermi regime (5 ^> 1, (3). Unlike in 
the present case, all vortices were found to be distinctly 
curved. Also, in Q the limiting {I = 1) frequencies for 
the axially symmetric vortex appear to coincide for both 
the S and the U-vortex. These differences should be at- 
tributed to the Thomas- Fermi nature of the condensates 
studied in the above references. 

In order to complete the study of S-vortices we note 
that the qualitative picture should be modified for low 



FIG. 6: Angular momentum as a function of the precession 
frequency for vortices in an axially symmetric elongated trap 
with p = 1/2, 5 = 0.28. The solid line shows the U-vortex 
branch. The dashed line is the S-vortex for which the £ = 
state (for uj < 0.26) is a dark soliton. Both modes become 
the axially symmetric vortex at the limit £ — > 1 but this is 
obtained at different frequencies. 



interaction strengths. This is because the solitonic vor- 
tex is a solution of the nonlinear theory which exists 
only above a certain value of the interaction strength, 
5 > Sq, which depends on (3 0, 0|. In order to ex- 
plore the regime 5 < 5a we choose 5 = 0.28 in the 
elongated trap with (3 = 1/2. Fig. shows the angu- 
lar momentum of the vortex modes as a function of the 
precession frequency. A substantial modification for the 
S-vortex is that now I — > for a finite limiting frequency 
uj = 0.26. For uj < 0.26 we obtain an axially symmetric 
wave function with a node on the z — plane which, 
in the context of the present nonlinear theory, is called 
a dark soliton by comparison to the dark soliton of the 
nonlinear Schrodinger equation. Summarizing, the pic- 
ture in Fig.Elcould also be described as following: a dark 
soliton which may initially be created in a trap would not 
respond to rotations of the system until we exceed a lim- 
iting frequency 0.26. For uj > 0.26 it could turn into 
an S-vortex which would become an axially symmetric 
one for uj > 0.49 (these frequencies should be compared 
with Eq. (0 in Sec. Ill when setting £ = and £ = 1 
respectively.) Finally, we note that the particle density 
plots for the S-vortices are qualitatively similar to those 
for higher interaction strengths presented above (as in 
Fig. [SJ, except, of course, for the £ = case. 

The U-vortex branch is qualitatively similar to its 
counterpart for 5 = 2 but its frequency range is shifted 
closer to the transverse trapping frequency (u>x) while it 
becomes narrower as 5 decreases. 

In Ref. ^3 a series of virial relations are derived which 
should be satisfied by any vortex solution. We have 
checked that the solutions presented here have, in partic- 
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ular, the counterintuitive property that the moments of 
the particle density vanish: J x \^\ 2 dV — j y\^\ 2 dV 
| z |^| 2 dV = 0. This is trivially satisfied by the S- vortex 
due to its symmetry along the z axis, but it is a nontrivial 
result for the U-vortex. 

We comment briefly on the relation of the vortex modes 
presented here to vortices in helium II. In Ref. an 
equation is derived for the dynamics of a vortex line in a 
rotating bucket containing the superfluid and it is found 
that there are vortex lines precessing on an unstable orbit 
around the bucket center at some distance r which is an 
increasing function of the precession frequency. Since the 
angular momentum of the vortex is a decreasing function 
of r (as discussed in 0), we conclude that the angular 
momentum is a decreasing function of the precession fre- 
quency, that is the situation is analogous to that for the 
U-vortex in Figs.EEl and© 

We finally turn to the delicate issue of stability of the 
calculated vortex states. According to 01 "mechanical 
stability" requires that the energy be a convex function 
of angular momentum: duo/dl = d 2 E/d£ 2 > 0. How- 
ever, this is a condition for stability of the solution in 
the rotating frame under variations that change the an- 
gular momentum. In experiments in which the trap is 
axially symmetric (i.e., when the system is no longer be- 
ing stirred), the angular momentum is conserved and this 
condition is irrelevant for the stability of the modes: what 
is relevant is that the mode should be dynamically sta- 
ble against decay into other modes which have the same 
angular momentum. As is evident from Fig. [21 the pre- 
cessing U-vortex has div/d£ < 0, so it does not satisfy 
the criterion of mechanical stability of Ref. [ljj. How- 
ever, the U-vortex must be dynamically stable since it 
is the lowest energy configuration for given angular mo- 
mentum < £ < 1 as we know from the calculations in 
the present section. Indeed the U-vortex has been seen 
in experiments in which angular momentum is con- 
served. 

It is also interesting to comment on the issue of overall 
thermodynamic stability of the rotating system. A single 
axially symmetric {£ = 1) vortex will have lower energy 
in the rotating frame than the nonrotating ground state 
if lo < Ey — E where Ey is the vortex energy and E 
is the ground state energy pj. Now, taking into account 
that lo = dE/d£, we may write 



E< 



En = 



u(t) dt. 



(9) 



Applying this relation for the U-vortex calculated for a 
spherical trap with <5 = 85 (see Fig.|2J we find Ey — E — 
0.23 or, in physical units, 1.8 Hz which lies in the range 
of precession frequencies for the U-vortex. 



III. DILUTE CONDENSATES 

We now turn to discuss more specifically dilute conden- 
sates which offer the possibility to study vortex modes in 



the limit of a small interaction strength where the energy 
levels of the system are close to those of the harmonic os- 
cillator. This analysis will provide an analytic des crip tion 
of the features of the S and U- vortices. In Ref. [l(| the 
energy E mt of a rotating system is minimised in the case 
of the very dilute limit in the sense that the nonlinear 
term in Eq. @ is much smaller than the other energy 
scales, namely ?iuj±, 7k<;||, or 



8< 1,13. 



(10) 



The problem is then reduced to the lowest Landau level 
(LLL) on the transverse plane (8 <§; 1) while the wave 
function in the longitudinal direction z is the Gaussian 
ground state of the ID harmonic oscillator. We refer to 
this as the 2D LLL regime. An analytic solution for vor- 
tices with angular momentum in the range < £ < 1 is 
given in |18|. The minimisation of E rot gives the Gaus- 
sian nonrotating ground state as the global minimum for 
lo < lojj = 1 — 8. For lo > lo\j the axially symmetric vor- 
tex becomes the global minimum of the functional JSJl. 
Exactly at ui = lojj there is a set of states with < £ < 1 
which are off-center vortices. These are denoted by a 
solid line in Fig. \7\ which shows the angular momentum 
as a function of the precession frequency. For higher val- 
ues of lo (but still lower than the transverse frequency of 
the potential) many-vortex states and vortex lattices be- 
come the absolute minimum 0, . All vortex lines are 
parallel to the symmetry axis of the trap since this is a 
2D calculation. For larger values of 8 these evolve in the 
curved vortex lines as, e.g., in Fig. ^ and thus they cor- 
respond to the U-vortex mode discussed in the previous 
section. 

In elongated traps excited states of the rotating con- 
densate can be found without mixing higher Landau lev- 
els. This is an interesting issue also because most ex- 
periments involving vortices, and in particular those in 
which the shapes of individual vortex lines were studied, 
were performed in very elongated traps, that is at (3 <C 1 

0. We first consider the case of a non- interacting gas, 

1. e., 8 = 0, whose first excited state has a node in the z 
direction: 



I") 



its energy is Ei = (3, and has vanishing angular momen- 
tum. We note that for increasing 8 the wave function 
evolves to the dark soliton in the framework of the non- 
linear Gross-Pitaevskii model. A further excited state is 
of the vortex type: 



*2 



1 



T 3/4 



(12) 



and it has energy Ei — 1 and angular momentum £ = 1 . 
It is thus obvious that for the frequency 



u s = 1-/3 



(13) 



the two states 111 111 and l|12|l have the same energy E TOt 
in the rotating frame. For lo < lo$ we have as the first 
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FIG. 7: The angular momentum i as a function of the fre- 
quency for i5< 1,(8,1-/3 for single vortex states. The solid 
line is the result of Ref. [Toll . The dashed line is a plot of 
Eq. I)lri|l . The horizontal axis is not to scale. The wave func- 
tion at I — for the dashed line (u) < 1 — j3 — 5) is while 
for the solid line (ui < 1 — S) it is the Gaussian. The 1 = 1 
state is ^2 for both cases. 



excited state while *f?2 becomes the first excited state for 
uj > u>s- 

The simple picture is modified in the presence of a 
weak interaction term, that is, for 5 small but finite. For 
6 <C /3, 1 — (3 we make the ansatz 



* = Ci*i +c 2 * 2 , 



(14) 



where the coefficients ci,C2 are constrained by \c\\ 2 + 
|c 2 | 2 = 1 and |c 2 | 2 = I. Substituting * in J3J we find 
the energy of the first excited state as a function of the 
angular momentum i: 

E(£)= 6 -l 2 + (l-P-5)£ + P+^ 1 (15) 

where we have subtracted the ground state energy of the 
harmonic oscillator 1 + (3/2, a convention that will be 
adopted from now on. Now, E mi {t) has a minimum for 
any I in the range < I < 1 and a frequency u> = dE/d£ 
or 



uj = 1 - (3 - 5(1 - £), 



(16) 



as depicted by the dashed line in Fig.[7| For uj < 1 — /3 — 5 
the wave function is ^ — ^1 (1 = 0) while for u> > 1 — /3 
we have * = * 2 {I = 1)- 

The vortex line represented by the ansatz l|14|) is found 
by setting 

* = => V2ci z + c 2 {x + iy) = 

y = 0, V2ci z + c 2 x = 0, (17) 

where for simplicity we may consider c\ , ci real. For 
frequencies in the range 1 — [3 < uj < 1 — [3 — <5 we 



have ci , C2 7^ and thus Eq. lfl7|l gives a straight vor- 
tex line which is tilted with respect to the symmetry 
axis z. Summarizing, the first excited state IjMfl has a 
wave function with a nodal plane perpendicular to the 
z axis for low frequencies which turns to a tilted vor- 
tex for 1 — (3 — S < uj < 1 — (3. This gradually aligns 
with the symmetry axis as the frequency increases. It 
becomes the axially symmetric vortex above the limiting 
frequency uj > uis — 1 — f3. This is the weak interaction 
limit of the S- vortex discussed in the previous section. 

Both the S and the U-vortex become the axially sym- 
metric vortex for a sufficiently high uj. However, the 
corresponding limiting frequencies for the two types of 
vortices do not coincide. In fact, there is no sense of pre- 
cession in an axially symmetric vortex and the limiting 
frequency may be absorbed into an effective chemical po- 
tential ji — fi + uj while the wave function reduces to a 
quasi static configuration with chemical potential ft. 

We are able to take our calculation one step further 
if we assume a very elongated trap, i.e., (3 <C 1, and 
an interaction strength comparable to the latter energy 
scale: 



8 < 1, and 5 - [3. 



(18) 



In this case it is useful to extend the LLL calculation to 
3D by including higher axial oscillator states. We use the 
wave function ansatz 

{pA)Uz), (19) 

where c m „ are complex constants, 



Xm(p, (f>) 



1 



1/2 



e' m *p m e- p!/2 (20) 



are eigenstates of the 2D harmonic oscillator and 



2" 



(2n- 1)! 



1/2 



fl/4 



z e 



-z 2 /2 



(21) 



form a basis for the expansion of the wave function in 
the axial direction. We aim to find the vortex states for 
a fixed angular momentum within the present ansatz. An 
efficient way to do this is by minimizing the Lyapunov 
functional 



^rot - E 



a 2 



b 2 f, 



(22) 



where I is the angular momentum functional of Eq. Q 
and v = J 1^1 2 dV . The constants ai , a 2 , 61, 62 are chosen 
so as to obtain a solution for the desired number of par- 
ticles and angular momentum. We substitute the wave 
function 1)19(1 in Q22JI and minimise E" ot with respect to 
Cmn by a conjugate gradient method. Computer limita- 
tions allow for modes typically up to m = 4 and n = 14. 

When the interaction term is small, i.e., 5 <C (3 the 
algorithm converges to the two single-vortex solutions 
which were discussed in Fig. [7| Specifically, for the U- 
vortex we have that all c mn with n ^ vanish. For the 
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FIG. 8: The angular momentum £ as a function of the fre- 
quency for the S (dashed line) and U- vortex (solid line) within 
the 3D LLL calculation for <5//3 = 2. Here, dl/dui > (except 
for small values of €) unlike in Fig. 



S- vortex we have that only coi,cio are nonzero. When 
S becomes comparable or greater than (3 several c mn are 
nonzero and their values decrease with respect to both 
indices. For the S- vortex we find that all c mn with m + n 
even vanish, and this is consistent with the symmetry of 
an S-shaped vortex. For the U- vortex we find that for all 
n odd c rnn — which is consistent with the symmetry of 
a U-shaped vortex. 

Fig.|S|shows the angular momentum of the two vortices 
as a function of precession frequency for parameter values 
8/(3 = 2. The most dramatic effect is seen in the U- vortex 
where d£/duj is now positive except for small values of £. 
This means that the abrupt transition from the £ — 
to the £ = 1 state depicted in Fig. [7| is smoothed out 
and the angular momentum of the U-vortex is now a 
smooth function of the frequency lo, although the range 
of frequencies corresponding to < £ < 1 is narrow. The 
vortices in the latter range are located off-center while the 
contribution of modes with n > in the wave function 
dl9fl means that vortex lines are curved. The S-vortex 
lines are also curved and the slope of the branch (dashed 
line in Fig. [SJ) is not constant unlike in Fig. [7\ 

A concise presentation of the main features of the two 
vortex modes can be given by the graph of the energy as a 
function of the angular momentum shown in Fig. EI The 
energy of the U-vortex mode for £ = is the energy of 
the nonrotating ground state. The corresponding value 
for the S-vortex is the energy of the wave function \I>i. 
The S-vortex has higher energy than the U-vortex but 
both curves converge to the same value for £ = 1 which 
is the energy of the axially symmetric vortex ^>2- 



FIG. 9: The energy E of Eq. @ as a function of the an- 
gular momentum £ for the S (dashed line) and the U-vortex 
(solid line) of Fig. GO where 5//3 = 2. The energies of the 
nonrotating ground state, and of the wave functions ^i, and 
\p2, are Eo, E\, and E2 correspondingly. (In the noninteract- 
ing system the corresponding energies are 0, (3, and I, having 
subtracted the ground state energy of the harmonic oscillator 
1 + /3/2.) 

IV. CONCLUSIONS 

We have analysed the two single vortex modes, the 
S-vortex and the U-vortex, which are stationary states 
of the Gross-Pitaevskii model precessing at a constant 
frequency u> in a BEC confined in an axially symmet- 
ric harmonic potential. Vortices for the whole range of 
angular momenta < £ < 1 have been studied. The 
precession frequency range is distinctly different for the 
two modes. The U-vortices fall within a relatively nar- 
row frequency range. For those U-vortices studied in Sec. 
fl we find that the precession frequency is a decreasing 
function of the angular momentum. The S-vortex is an 
excited state in elongated condensates. Their frequency 
range extends down to zero (for sufficiently high interac- 
tion strengths) and the nonrotating member of the fam- 
ily is the so-called solitonic vortex which is a solution 
of the nonlinear Gross-Pitaevskii model 0, . Both the 
U-vortex and the S-vortex have been directly seen in ex- 
periment 0, Q ■ Perhaps the S-vortex is also related to 
shape deformations of vortex lines observed in the form 
on Kelvin modes 

A complete analysis for vortex solutions could be car- 
ried out for a dilute condensate within a 3D lowest Lan- 
dau level calculation which includes the axial oscillator 
states. We have given a derivation for the S-vortex which 
is complementary to the LLL calculation for the U-vortex 
10}. fncluding the axial oscillator states causes the tran- 
sition from a nonrotating ground state to an axially sym- 
metric vortex to happen smoothly as the rotation fre- 
quency increases: vortices with angular momentum £ less 



8 



than unity are thermodynamically stable. 

Both the S and the U-vortex become axially symmet- 
ric in the limit I — > 1 but this happens at different fre- 
quencies in the two cases. This seems to be in disagree- 
ment with the results reported in |7[, however, one should 
note that in the latter reference very elongated and very 
dense condensates are studied. Although the results are 
related, an extrapolation of our results to the Thomas- 
Fermi regime, and thus a direct comparison to Q, does 
not appear straightforward. We further find that, for 
weak interactions, an S-vortex line is only slightly bent 
in contrast to the original observations in [2J. In fact, 
both vortex modes in the dilute limit presented in Fig. [3 
represent straight vortex lines. However, as the density 
increases vortex lines tend to bend (as, e.g., in Figs. l4l5jl . 



Thus, the difference between the present results and the 
observations is clearly due to the difference in the density 
of the condensates. 
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